hep-ph/ 0306071 
IFIC/03-24 
ZU-TH 07/03 
UWThPh-2003-05 

Testing the Mechanism of R-parity Breaking with 
Slepton LSP Decays 



A. Bartl 1 , M. Hirsch 2 , T. Kernreiter 1 - 2 , W. Porod 3 and J. W. F. Valle 2 

1 Institut fur Theoretische Physik, Universitdt Wien, 
A- 1090 Vienna, Austria 
2 AHEP Group, Instituto de Fisica Corpuscular - C.S.I.C. 
Universitat de Valencia 
Edificio de Institutos de Paterna, Apartado 22085, E-^6071 Valencia, Spain 
3 Institut fur Theoretische Physik, Universitdt Zurich, 
CH-8057 Zurich, Switzerland 



PACS: 12.60 Jv, U.60Pq, 23.40-s 

Abstract 

In supersymmetric models R-parity can be violated through either 
bilinear or trilinear terms in the superpotential, or both. If charged 
scalar leptons are the lightest supersymmetric particles, their decay 
properties can be used to obtain information about the relative im- 
portance of these couplings. We show that in some specific scenarios 
it is even possible to decide whether bilinear or trilinear terms give 
the dominant contribution to the neutrino mass matrix. 



1 Introduction 



The most general form for the R-parity violating and lepton number violating 
part of the superpotential is given by 



Taking into account the asymmetry of Xijk in % and j, Eq. contains nine 
different Ay&, 27 A^- fc and three e,, for a total of 39 parameters. In principle 
- as already pointed out in 1J - the number of parameters in Eq. can 
be reduced by 3 by a suitable rotation of basis and a number of authors 
have used this freedom to absorb the bilinear parameters q into re-defined 
trilinear parameters [2J. 

However, for consistency one has to add three more bilinear terms to 
V so ft, the SUSY breaking potential, 1 



It is important to note that there is no rotation which can eliminate the 
bilinear terms in Eq. (j2J) and Eq. simultaneously 1J 2 and thus we are 
back to 39 independent parameters. 

With this huge number of parameters it seems hardly possible to make 
any definite prediction for phenomenology. Could one reduce the number 
of free parameters in the theory? On theoretical grounds one could argue 
that a model with only bilinear terms is a self-consistent theory, whereas a 
model with only trilinear terms is not. The reasoning behind this argument 
is quite simple: Renormalization group running |Hj will generate bilinears 
whenever trilinears are present, but RGEs which start with bilinears only 

1 We will not consider the corresponding lepton number violating trilinear soft masses, 
Aijk and A' i j k , since they will not affect our conclusions. 

2 Unless the parameters Bi and the corrsponding MSSM parameter B are equal (as well 
as m 2 Li = rn 2 Hd ). However, such an equality is unstable under RGE running [3]. 
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will never generate trilinears [3]. In addition, models which break R-parity 
spontaneously [5] at the TeV scale [H] will lead to effective low-energy models 
where only bilinear terms are present. 

On the other hand, experiments are done in the physical mass-eigenstate 
basis. Before exploring the implications of Eq. (JTJ, one therefore has to 
re-diagonalize all mass matrices of the model. The essential point here is 
that the bilinears introduce mixing among the various states of the model 
and thus lead to the appearance of "effective" trilinears in the mass basis, 
for example, 

K33 ~ — h b (3) 
A 4 

where hf, and \x are the bottom Yukawa coupling and the Higgs mixing param- 
eter. These "effective" trilinears are unavoidable even in models which start 
with bilinear terms in the superpotential only. The essential point to realize 
here, however, is that these "effective" trilinears necessarily always follow 
the hierarchy implied by the Yukawa couplings of the standard model and 
therefore are not new free parameters of the theory in contrast to "genuine" 
trilinears. 

All terms in Eq. (0) violate lepton number by one unit and thus will 
necessarily contribute to the (Majorana) neutrino mass matrix. Quite a 
number of articles have investigated the consequences of Eq. ([TJ for neutrino 
physics, some of them considering only the bilinear terms [Zj , [Hj others only 
the trilinear terms and a few have also entertained the possibility that 
both kind of terms exist [10 . 

An obvious question to ask then is: Which of the terms in Eq. ( f l j) 
give the dominant contribution to the neutrino masses? And, is there any 
possibility to settle this question experimentally! The aim of this paper is to 
demonstrate that, if charged scalar leptons are the lightest supersymmetric 
particles (LSPs), there are different broken R parity scenarios where one can 
probe for the origin of neutrino mass. 

How likely is it that charged sleptons are the LSPs? In models of su- 
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persymmetry breaking based on minimal supergravity boundary conditions 
(mSugra) jTT] , in gauge mediated SUSY breaking jT2j and in anomaly medi- 
ated SUSY breaking [T^ITl] this actually happens in sizeable portions of the 
parameter space. Usually this possibility is declared "ruled out cosmologi- 
cally" and simply discarded (except in GMSB, where the charged sleptons 
are actually the NLSPs). However, with R-parity violated the charged scalars 
decay, so cosmological limits simply do not apply 3 

It is exactly the observation that the LSP decays through the same Im- 
parity violating operators, which also govern the entries in the neutrino mass 
matrix, which forms the basic idea of the current paper. To give a trivial 
example, in a world were only A^'s are non-zero, charged scalars can decay 
only to leptonic final states, whereas if only A^- fc 's are non-zero, there will 
only be hadronic final states. Of course, this example is vastly oversimplied, 
but here we will argue that, with experimental information from neutrino 
physics [15 , it is possible to fix the size of the various couplings assuming 
one or the other combination being dominant and check the consequences 
for the charged slepton decays. 

Upper limits on all parameters in Eq. (jTJ exist and have been extensively 
discussed in the literature 0EE], jXZOHl - However, in agreement with [12], we 
find that limits on the couplings imposed by current neutrino oscillation data 
are usually much stronger than all other indirect limits. 4 Notable exceptions 
from this general rule are the limits from double beta decay [IBj, from non- 
observation of n — > 3e, fiTi — > eTi, and from Attlb and ArriK Here 
we note only that none of the existing limits is strong enough to require 
charged scalar leptons decaying with visible decay lengths. 

Scalar tau LSPs have been discussed in the literature before. In |2Hj scalar 
tau decays within bilinear R-parity violation have been discussed for a one 
3 A possible dark matter candidate in this case is the axion. 

4 This claim is true only if the left-right mixing in the squark/slepton sector is not 
exactly zero. Exactly vanishing left-right mixing can not be excluded at present, but 
requires fine tuning of parameters. We will not consider this possibility. 
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generation model of R-parity violation. In [21] possible similarities between 
the phenomenology of a charged Higgs and a stau have been discussed within 
trilinear R-parity violation. In [22J charged scalar LSP decays within bilinear 
R-parity violation (with 3 generations) have been discussed with special em- 
phasis on their relation to neutrino physics. 5 And in GMSB the scalar tau 
can be a quite long-lived next-to-lightest supersymmetric particle (NLSP), 
which from the collider point of view looks like a stable, charged LSP, see 
for example ^2]- Finally charged scalar LSPs have been discussed in an R- 
parity conserving extension of the MSSM in [21]. These authors argue that 
one should not take into account the cosmological limits on charged LSPs, 
because alternative cosmologies can provide loopholes to the usual "overclo- 
sure" argument j'JI) . 

This paper is organized as follows: In the next section we will set up the 
model definitions. In section 3 the charged scalar decay widths are given. 
Section 4 summarizes some necessary formulas for the neutrino mass ma- 
trix, while section 5 discusses charged scalar lepton decays under different 
hypotheses both analytically and numerically. We then close with a short 
summary and outlook. 

2 The Model 

We work in the minimally supersymmetrized version of the standard model 
(SM) [H], augmented with the R-parity violating terms in Eq. Super- 
potential and mass matrices of the MSSM have been given exhaustively in 
the literature and will not be repeated here. For mass matrices including the 
bilinear R-parity violating terms see, for example [Zj. 

As will be seen in the next section, the slepton decay widths depend on 

5 After completion of this article, the preprint |23| appeared. The authors study general 
R-parity violation in mSugra considering also the case of a scalar tau LSP. However, the 
emphasis of |23| is on the supersymmetry breaking parameters and no attempt is made to 
explain current neutrino data by the R-parity violating parameters. 
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left-right mixing in the scalar sector. We therefore give a short repetition of 
£i - £r mixing. The masses and couplings of the £ follow from the symmetric 
2x2 mass matrix which in the basis (£r,£l) reads [TT|l2"H] 



/ M? M~ 

4, = -(4 A) lRR lRL ] ( : ] , (4) 

where 



)l 









Mj LL = M? +m| cos 2(3 ^sin 2 9^ - ^ ) + m 2 £ , (5) 



Mj RR = Ml -m 2 z sm 2 e w cos2(3 + m 2 i , (6) 
M L = M l =m t (A e - fittm P). (7) 



RL ZLR 



Here tan/3 = v u /vd with Vd(v u ) being the vacuum expectation value of the 
Higgs field H%(H®), ran is the mass of the appropriate lepton and Qw is the 
weak mixing angle, [i is the Higgs-Higgsino mass parameter and Mi, M^, An 
are the soft SUSY-breaking parameters of the slepton system ^T] . The mass 
eigenstates £i are (ii,^) — (^Ri^l)^ with 



cos 9 c sin Q- 



sin 9g cos 9~ t 



with 

M| Ml - m 2 

cos 9-, = , sin 9 



(Ml J 2 + (ml - MU J(Miy + K - M lf 

(9) 



ZLR' «1 t-RR' U l LR «1 «Hfl 
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while the mass eigenvalues are given by 

m t = \ (<*L + M L ) =f \/( M L- M U 2 + « M i/) • < 10 > 

Due to the appearance of the lepton mass in Eq. (J7|) one expects the left-right 
mixing is tiny in the seletron, small in the smuon and potentially sizeable in 
the stau sectors. 

In the R-parity violating version of the MSSM a priori any SUSY particle 
can be the LSP. However, supplementing the MSSM with mSugra motivated 
boundary conditions, one usually finds that either the lightest neutralino or 
one of the charged scalars is the LSP (CSLSP). As a rule of thumb, charged 
scalars are LSPs if m <C Mu 2 and fi large. CSLSPs in these models are 
usually mainly right sleptons. Furthermore, even though one expects some 
splitting between f , \x and e from RGE running, the latter are not much 
heavier than the former and so might also dominantly decay through R- 
parity violating operators even though they are not strictly the LSP. 6 

3 Decay Widths 

Charged scalar leptons lighter than all other supersymmetric particles will 
decay through Jji v vertices. Possible final states are either l{i>j or qq'. We will 
present the decay widths taking into account only the trilinear parameters. 
For the corresponding formulas for the bilinear terms see |22j . 

The relevant Lagrangian to study these two-body decays of l\ is obtained 
from Eq. (JTJ). It is given by [2|, 

Ant = e j ie k [(sm9e j Xijk)PL + (cose^Xik^PR^i 

~~ V *Ak Sil1 Qe&ldkPbUn + h.C. , (11) 
3 

6 This is similar to what happens quite generically in GMSB models, where such a 
scenario is called co-NLSP |12| . 
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where Pr,l = (1 ± 7s)/2 and V n j is the corresponding element of the CKM 
matrix. Here we work in a basis where the CKM matrix is solely due to the 
mixing of the up-type quarks, [efj-i denotes the lighter slepton mass eigenstate 
with family index j\. In the limit <C the leptonic two-body decay 
widths of £i read, 

T(eji -» e fc ^z/j) = ^[(sin O^Xjjk) 2 + (cos g . A ifci ) 2 ] (12) 

i i 

For £1 ~ £r, we find for the branching ratios 

5r (l,2,3) " 2 

where we have introduced the following shorthand notation: Br^ 2 3 ^ = 
Br(£i — > (e, /i, r) z/j). Corrections to Eq. (fT^j) are oc #| 1. In the 
Appendix we give the solutions to Eq. (fT^J) with respect to the trilinear cou- 
plings A 2 j7c by making an expansion in the small parameters 0$. In the limit 
01 — > , we derive form Eq. (|T5j) that 

Srf < 0.5, V i (14) 

It is interesting to note that Eq. (JTljl is a definite 'prediction (correct up to 
order #|) of trilinear R-parity violation despite the large number of param- 
eters. An experimental result contradicting Eq. (|14jl would be a clear sign 
for bilinear R-parity breaking, if the CSLSPs are mainly right (isosinglet) 
sleptons. 

In order to determine A 2 A a measurement of the total width of the slep- 
tons is necessary. Information on the widths could be obtained if a finite 
decay length is observed, we will return to this point later. However, be- 
cause the next to leading term in Eqs ([34 p -([36| ) is only of the order #| ratios 
A? fc /AL/ & may be determined without this information from the branching 



1 - 



(\2 \2 \2 \ 

i 23i' i3i' \jMl 
V A 2 



(13) 
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ratios Br{e k — > q ^2 i/ m ) alone. 

Note also that there are 9 different \j k and 9 observables, which implies 
that in principle it is possible to determine all A^ if there is no bilinear 
R-parity breaking. Conversely, adding also the bilinear terms in the super- 
potential and/or in the soft SUSY breaking potential a full reconstruction of 
all parameters from the leptonic decays of the CSLSPs is impossible in princi- 
ple. Differentiating between bilinear and trilinear R-parity breaking therefore 
necessarily requires the construction of conditions such as Eq. fll4j) . 

For the hadronic two-body decay widths of i\ one finds 



(rl) nk = rd -> u n d k ) = n c (3^- sin 2 e g | ^ v nj \' ijk \ 2 , (15) 



where N c = 3 is the number of colours. (3 = 1 for j ^ 3 and (3 = (1 — (^) 2 ) 2 
for j = 3. For the expected small mixing in the first generation of sleptons 
the hadronic width is highly suppressed compared to the leptonic width if 
A' < A. 

4 Neutrino Mass Matrix 

The Majorana mass term for the neutrinos is of the form 



where Vu denote the left-handed weak eigenstates and M v is a complex 
symmetric 3x3 matrix |2fij . 

Contributions to the neutrino mass matrix are induced at tree-level by 
the bilinear lfi p terms. They can be expressed as |2*7ll2^] 



3 



c 



mass 



- (v Li ) T C (M„)ui v Li > + h.c. 



(16) 



M l9 2 + M 2 g' 
4 det(A^ x o) 



/ A 2 A e A M A e A T \ 
A e A M A 2 A M A r 
\A e A T A M A T A 2 J 



(17) 
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where Aj = \xvi + v^i- Due to the projective nature of the effective neutrino 
mass matrix m e g , only one neutrino acquires mass at the tree level. Therefore 
m e g is diagonalized with only two mixing angles which can be expressed in 
terms of Aj: 

tanfl 13 = - Ae , (18) 

(A2 +A2)2 

tan£ 23 = -^. (19) 

One-loop contributions to the neutrino mass matrix from bilinear terms have 
been discussed extensively, see for example j7| and will not be repeated here. 
At 1-loop level there are also contributions due to trilinear ftp couplings. The 
relevant A-terms are given in Eq. (fTTj) . for \'^ k type couplings we need: 

C y D X' ijk [d jL d k P L Vi + d\ R vtP L d 3 \ + h.c. . (20) 

The full mass matrix from trilinear terms is then given by m 1 ~ loop = m x +m x ' 



where 



A(A') _ l(iVc) A (') , (') 
"Hi' - 327] -2 A ijk A i'kj 



( m 2j \ ( m 'lk 

— 7T- + rrij sin 29). In ' 
ml J 



"XjJ \ m lk 



(21) 

where is the appropriate fermion mass, 9j denotes the appropriate sfermion 
mixing angle and m\ 2 ^ are the corresponding sfermion masses. Note the 
manifest symmetry of this matrix as well as the presence of logarithmic fac- 
tors. 

We stress that Eq. (}2"T| is an approximation to the full 1-loop calculation. 
Numerically we have found that Eq. (|21j) is good only up to factors of ~ 2 
in extreme cases. In the numerical part of this work we therefore always 
diagonalize the full 1-loop corrected (7 x 7) neutrino-neutralino mass matrix. 
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5 Neutrino Physics versus Collider Physics 



In the following we consider some specific scenarios with bilinear and trilinear 
R-parity breaking terms being present in different combinations and discuss 
the resulting decay patterns of the charged sleptons. We will calculate the 
decay lengths 



assuming a centre-of-mass energy y/s = 0.8 TeV, as typically could be 
reached in a future linear collider [SHUSH!- 

In order to reduce the number of parameters, the numerical calculations 
were performed in a constrained version of the MSSM. We have scanned the 
parameters in the following ranges: My 2 G [0,1.2] TeV, G [0,2.5] TeV, 
m G [0,1.2] TeV, A /m and B /m G [-3,3] and tan/3 G [2.5,10]. All 
randomly generated points were subsequently tested for consistency with the 
stationary conditions of the scalar potencial as well as for phenomenological 
constraints from supersymmetric particle searches [SU- To be conservative 
we require the charged scalars to be heavier than 100 GeV, although existing 
limits are somewhat weaker and depend on the flavour of the charged scalar 
|31j . We then select points in which charged scalars are the LSPs. This latter 
cut strongly prefers m <C M x / 2 and large values of /i. Note again that our 
charged scalars are mainly right sleptons, as discussed above. Also the low 
values of tan f3 in our scan lead to masses for the three different generations of 
charged sleptons which are rather similar such that all three charged sleptons 
decay dominantly through R-parity violating decay modes. 

The ftp parameters were chosen in order to fulfill the requirements of the 
various scenarios to be discussed and as input from neutrino physics we have 
used, unless noted otherwise, the following ranges, as currently determined 
from neutrino oscillation experiments [T5] : 




(22) 



0.3 < sin 6j\ tm < 0.7 



(23) 



11 



1.2 x 1(T 3 eV 2 < Ami 



Atm 



< 4.8 x 10" 3 eV 2 



(24) 



0.29 < tan 2 m < 0.86 



(25) 



5.1 x 10~ 5 eV 2 < Am 2 < 1.9 x 10" 4 eV 2 



(26) 



sin 2 13 < 0.05 



(27) 



Before starting the discussion, let us briefly summarize the extreme case 
in which neutrino oscillation parameters are determined only by bilinear pa- 
rameters [22]- Full details of this particular situation can be found in [22] 
and will not be repeated here. We will simply mention the following features 
which are the most important tools to distinguish the pure bilinear case from 
more complicated scenarios: 

• The hierarchy in Yukawa couplings implies that the slepton decay 
lengths must be very different, L(fi) <C L(fli) <C L(e\), in particu- 



• Hadronic final states have branching ratios which are much too small 
to be observable 

• The ei decays dominantly to the final state e v (electrons plus miss- 
ing (tranvserse) momentum) 

Next we will discuss two scenarios with non-zero trilinears. 



lar L(fx)/L(/ii) ~ (^) 2 



Ratios of branching ratios must obey the following condition: 
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5.1 Scenario I: \' ijk < X ijk 

Assume the atmospheric mass scale is generated by m e g in Eq. (j!7)) , whereas 
the solar mass scale is due to m A in Eq. ()21j) . If the couplings Ayi and Xij2 
are not much larger than Ay3 7 the leading contribution to m^, comes from 
the fj — r loop and the contributions of all other l± — t loops are subleading. 
Then m\, is given approximately as: 



, 1 / m| 

m 7H -m T sin 26 f In — =- 



/ A^ 33 A133A233 o\ 

A133A233 A233 

y oj 



(29) 



If the parameters in Eq. (|29|) are indeed such, that the solar mass squared 
difference is correctly explained, the scalar tau will decay with a very short 
decay length, as demonstrated in Fig. ^ Note that for a possible linear 
collider currently one expects to be able to measure decay lengths down to 10 
fim (32J in an event-by-event analysis. For typical expected luminosities [33] 
and production cross sections for charged scalars as calculated in [22] the 
expected number of events is about 10 4 events per year, implying that around 
10 events decay with length ~ 10 /im if its mean value is ~ 1.5 //m. In view 
of the small decay lengths implied by this scenario (Fig. we think that 
attempts by experimentalists at improving the measurement of very short 
decay lengths would be very interesting. 

For other A^-fc to provide the solar mass scale, the decay lengths of the 
corresponding sleptons would have to be even shorter. In this sense, scenario 
I is a "worst-case" scenario, since only ratios of A^-fc - and a lower limit on 
the absolute values of these parameters - could be measured. However, it is 
possible to turn this argument around and state that, if a finite decay length 
for ei and fti is found Ajji and A^ can not contribute significantly to the 
neutrino mass matrix. 



7 Whether this is indeed the case may be checked, in principle, by a comparison of the 
decay length of the sleptons. 
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I I 

10- 10 10" 9 10" 8 10" 7 10" 6 10" 5 10" 4 

Am% [eV 2 ] 

Figure 1: Stau decay length times ( j^^y ) 2 m [^m] versus solar mass squared 
difference Am^ [eV 2 ] in scenario I, defined as Xijk being responsible for the 
solar mass scale. If the trilinear loop is responsible for the solar mass, the 
stau will decay with a length too small to be measured. 

With a decay length of f\ too small to be measured, only a consistency 
check of this scenario is possible. This check is provided by the measure- 
ment of a certain ratio of branching ratios, which is correlated with the solar 
angle, as demonstrated in Fig. (J2J). We want to stress, however, that this 
measurement can not distinguish trilinear terms from a bilinear-only world, 
since the latter leads to a very similar prediction for the branching ratios. 
Instead only disagreement between prediction and measurement would es- 
tablish that neither scenario I nor bilinear loops are responsible for the solar 
neutrino mass scale. 

5.2 Scenario II: \jk <C A^- fe 

We now turn to the case where X^k <C A- jfc . If the A- jfc do not follow an inverse 
hierarchy, then the main contribution of the trilinears to the neutrino mass, 
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Figure 2: Ratio of branching ratios (1 — 2 x Brl^fil — 2 x Br"[) versus the 
solar angle, tan 2 9 & . The figure to the left is for VA^ non-zero and of similar 
magnitude, while the figure to the right assumes A 123 = 0. 

Eq. (jUJ), is approximately given by 



m « 



16n 



3 / m i 

nib sin 29~ h In ' 



1 A 133 
^133^233 



\^'l33^3 



333 



\' V \' V \ 

^133^233 A 133 A 333 > 



X' 2 

A 233 

^233^333 



y y 

/V 233 /V 333 



A' 2 



333 / 



(30) 



since the hi — b loop is largest due to the larger left-right mixing expected in 
the sbottom sector and rrid <C m s <C m&. Note that m x ' has also a projective 
structure, similar to m eff in Eq. (fT7j). but with different parameters. 
There are two limiting case, which we will study in some detail: 

• Scenario IIA: Bilinear R-parity breaking responsible for solar physics, 
A^ 33 for the atmospheric mass scale 

• Scenario IIB: Bilinear R-parity breaking responsible for atmospheric 
physics, A^ 33 for the solar mass scale 

Consider scenario IIA first. Here, due to the special structure of Eq. (}30|) . 
tan #23 and tan #13 can be approximately expressed as 



tan 9- 



23 



A^ 

a; 



233 



(31) 



3.33 
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IAWA333I (11)33/01)33 x 10 3 

Figure 3: Ratio A233/A333 (left panel) and (1^)33/(1^)33 (right panel) versus 
the atmospheric angle, tan 2 0Atm, for data points satisfying the criteria de- 
fined as scenario IIA. The plot to the right contains only points in which the 
decay ib is kinematically possible. 



tan 13 « A/l33 . (32) 

V^) 2 + (A 333 ) 2 



The ratio in Eq. (}3*Tj) squared can be related to the observable 

01)33 „ A 233 y w m; ^ 



(33) 



where (1^)33 denotes the partial width of t\ decaying into the final state ib. 
Due to the tiny selectron mixing we expect the number of events N(ei — > 
ib) <C N(p,\ — > ib), such that the corresponding relation for Eq. (}32|) will be 
impossible to test. 

For simplicity let us first consider the case in which only the A^ 33 are 
non-zero. Modifications of the results for other \'^ k 7^ will be discussed at 
the end of this section. 

In Fig. El we show the atmospheric angle versus A233/A333 (on the left 
panel) and (T^) 33/(^1)33 (on the right panel) demonstrating to which accu- 
racy one expects Eqs. (}3~T]) and (|3^|) to work. Obviously Eq. (J3~3*|) can be 
tested only if is bigger than m t + m^. 

Fig. IH (to the left) then shows the calculated decay lengths for fi, /ti 
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Figure 4: Decay lengths for fx, fix and ex, L(fx), L(flx), L{e~x) 
[fj,m], (bottom to top) versus the mass of the decaying particle, 
mr [GeV]. On the left panel: Scenario IIA, right panel scenario IIB. The 
plots assume that only X' i33 are different from zero, see text. 



and ei in scenario IIA. The top threshold for the f\ and fix decays is clearly 
visible, while it is invisible for the selectron case due to the negligible left- 
right mixing dictated by the small electron Yukawa coupling. This explains 
why the e\ decay length does not vary much over the range of masses shown. 
This is due to the fact that the decay of the e~x is completely dominated by 
the final state e^z/, induced by the bilinear term |A|. For the f\ below 
threshold the decay is dominated by the final state cb, which is induced from 
A 333 due to the non-zero value of V c f,, with a branching ratio into being 
a few percent. Above threshold the branching into ib quickly grows to more 
than 99 %, with the final state cb being suppressed by a factor of V^ b and the 
final state r v is tiny. For the decay of the fix, on the other hand, below 
threshold /i v dominates, with a branching ratio into cb of about ~ (1-20) 
%, while above threshold tb becomes quickly dominant, but (some) percent 
of the decays still go to /i v. Decays en — > ^2 v for % 7^ j are always very 
small, due to the assumed smallness of Ay&. 

From Fig. H] one concludes that an experimental test of Eq. ()33|) might 
be possible, but the decay lengths of the f\ are so short such that either 
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(i) an improvement in detecting short decay lengths will be necessary or (ii) 
additional input from the sbottom sector is needed. 8 

Turning to scenario IIB, in the right panel of Fig. E]we plot the calculated 
decay lengths of the charged scalar leptons as a function of the charged scalar 
mass. The main difference to scenario IIA is that the f\ now decays with a 
visible length, apart from some exceptional points where is very close to 
\/s/2. This can be traced to the fact that smaller values of A -33 are needed in 
this scenario since the solar mass scale is smaller than the atmospheric one. 
Note that the plot assumes that all A^ 33 are of the same order of magnitude. 
While A333 can not be much larger than the values used in this plot without 
violating the main assumption of scenario IIB, it could well be smaller than 
A' 133 or A233, since non-zero values for the latter are sufficient to create a non- 
zero solar neutrino mass difference (and mixing). We note in passing that in 
the (unlikely) case of X' 333 = the fi decay lengths would approach the e~\ 
decay lengths shown, with the dominating final state being (again due 

to the assumed absence of A^-3) . 

Let us finally discuss modifications of the results discussed for scenarios 

IIA and IIB once we allow for other \\- k 7^ 0. Non-zero X'^ k would, of 

course, affect the branching ratios discussed above. Since final states 

can not be distinguished experimentally for j = 1 and k = 1,2, only the 

sum of the corresponding (A^ fe ) 2 are measurable. More important for us, 

however, is that additional final states will shorten the decay lengths of the 

sleptons compared to the ones shown in Fig. 0] Numerically we find that 

in scenario IIA the decay length of the /2i is always visible unless X' 2 j k for 

j, k 7^ 3 are several times larger than the values for A 233 needed to explain 

the atmospheric neutrino mass difference. The situation is different for f\. 

Here, values of A 3jfc for j,k ^ 3 smaller than the "correct" A 333 by a factor of 

a few are needed, otherwise the decay length of the fi becomes to short to be 

8 Fixing the sbottom mixing angle and masses would allow to calculate the absolute size 
of A333 needed for scenario IIA, thus fixing the total fi width. 
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realistically measurable. In scenario IIB similar comments apply, although 
in this case even larger ratios AL-^/A^ are allowed for a visible decay length 
of the \x\ decays. For the f\ in scenario IIB A 3A . as big as A 333 are possible 
for decay lengths to be observable. The decay of the ex, on the other hand, 
never shows any sensitivity to A' ljfc , again due to the tiny left-right mixing in 
the selectron sector. 

6 Conclusions 

We have studied the decay properties of a right charged scalar lepton LSP 
(CSLSP) in a general model of R-parity violation containing both bilinear 
and trilinear terms. Branching ratios and decay lengths of CSLSPs contain 
information about ratios and absolute values of R-parity violating couplings 
which at the same time contribute to the neutrino mass matrix. We have 
investigated to what extent it is possible to test experimentally whether 
bilinear or trilinear terms give the dominant contribution to the neutrino 
masses. 

Due to the huge number of parameters characterizing trilinear models 
the resulting phenomenology can be quite diverse. In fact one can envis- 
age a mixed situation in which bilinear terms are responsible for generating 
the neutrino masses required to account for current oscillation data, yet the 
CSLSP decays are governed by trilinears. Thus in general the complexity 
of the physics resulting from Eq. (jlj is such that the existence of non-zero 
bilinear terms can not be established. 

Nevertheless, we have shown that different scenarios for neutrino masses 
exist which lead to considerably different phenomenology and thus allows 
for the origin of neutrino mass to be experimentally probed. In some cases 
we find that bilinears and trilinears can be clearly distinguished. Especially 
noteworthy is that in trilinear-only models all right CSLSPs should obey 
Br(£i — > (e, /x, r) ^ z/j) < 0.5. This is to be contrasted with the bilinear 
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model predicting [22] Br{e\ — > ej^z/j) ~ 1. 

Decay lengths for the CSLSPs could be measurably large or very small 
depending mainly on the absolute values of A^. Observing a finite length 
for §i and/or fix would establish that the corresponding A^i and Xij2 do not 
contribute significantly to the neutrino mass matrix. If Ay& are somewhat 
smaller than \[- k hadronic final states will have measurable branching ratios, 
at least for the f\ and depending on X^/X^jk also for the fx. In contrast, note 
that in the bilinear-only model |22] hadronic final states are never visible. 

Finally, since reasonable scenarios for the neutrino mass matrix exist in 
which the decay length of the fx is just at or below the borderline of what 
is currently thought of being experimentally accessible [30J we stress that 
efforts to optimize decay length measurements might be worth undertaking. 
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7 Appendix 



Here we give the solutions for the trilinear couplings Af jfe by inverting Eq. f!12[) 
and making an expansion in 9$ up to 0(9~). This yields 

A? 2 i = ^ C- 1 {rf ot (l - 25rf) + 0l(rf ot (5rJ - Srf) - lf ot £rf + Tt ot Brf)}, 
A?3i = 5 C*" 1 {rf ot (l - 25rf) + 9l(Tt ot (Brl - Brl) + I&flr? - T^Brf)}, 
A 2 23 i = J C- 1 {rf ot (l - 2£rf) + ^(rf ot (Sr| + 5r 2 g ) - I&flrf - 17^)}, (34) 



2 
1 

2 



AL = i C" 1 {rf ot (l - 2Br{) + ^(ittCSrf - Brl) + T^Br* - r[ ot 5r 2 f )},(35) 



A? 23 = \ C- 1 {TIM ~ 2Brl) + 9 2 M ot {Br\ + Brl) - ^4 - T^Brj )}, 
A? 33 = ^ C" 1 {ltot(l " 2Sr|) + ^(r[ ot (5r[ - Br|) - rf ot £rf + I^Br? )}, 
AL = J C- 1 {!**(! - 25r[) + 9l{Y\ ot {Brl - Brl) + ^4 - I&Srf )},(36) 



where r^ ot denotes the total leptonic decay width of the appropriate scalar 
lepton, which is, for the case where l\ ~ £r, not much different from the 
total decay width if A' < A. C = where we have made use of the fact 
that — rrifa ~ m fl (= m^). As can be seen in Eqs. (J54*j) - (|3l)j) the next to 
leading term is only of the order 9 2 . The measurement of I^ ot and 9~ t provides 
the information whether the expansion to 0(9j) is sufficient. 



21 



References 



[1] L. J. Hall and M. Suzuki, Nucl. Phys. B 231 (1984) 419. 

[2] For a review on trilinear R-parity violation, see: H. K. Dreiner, 
|hep-ph/9707435| 

[3] B. de Carlos and P. L. White, Phys. Rev. D55, 4222 (1997), 
lhep-ph/9609"443| |; B. de Carlos and P. L. White, Phys. Rev. D54, 
3427 (1996), |hep-ph/960238l"| ; E. Nardi, Phys. Rev. D55, 5772 (1997), 



lhep-ph/9610540| . 

[4] M. A. Diaz, J. Ferrandis, J. C. Romao and J. W. F. Valle, Nucl. Phys. 
B590, 3 (2000), lhep-ph/9906343 . 



[5] C. S. Aulakh and R. N. Mohapatra, Phys. Lett. B119, 136 (1982). G. G. 
Ross and J. W. F. Valle, Phys. Lett. B151, 375 (1985). J. R. Ellis et 
al, Phys. Lett. B150, 142 (1985). 

[6] A. Masiero and J. W. F. Valle, Phys. Lett. B251, 273 (1990). J. C. 
Romao, C. A. Santos and J. W. F. Valle, Phys. Lett. B288, 311 (1992). 
J. C. Romao, A. Ioannisian and J. W. F. Valle, Phys. Rev. D55, 427 
(1997), |hep-ph/960740T| . 

[7] M. Hirsch, M. A. Diaz, W. Porod, J. C. Romao and J. W. F. Valle, 
Phys. Rev. D62, 113008 (2000), |hep-ph/0004Tl5] ; [Erratum-ibid. D 65 
(2002) 119901]; M. A. Diaz, M. Hirsch, W. Porod, J. C. Romao and 
J. W. F. Valle, D68, 013009 (2003) |hep-ph/030202l| 

[8] R. Hempfling, Nucl. Phys. B478, 3 (1996), hep-ph/ 95 112881 ; D. E. 
Kaplan and A. E. Nelson, JHEP 01, 033 (2000), lhep-ph/9901254| ; S. Y. 
Choi, E. J. Chun, S. K. Kang and J. S. Lee, Phys. Rev. D60, 075002 



(1999), 


hep-ph/9903465 


; E. J. Chun and S. K 


075012 (2000), 


|hep-ph/9909429 


; E. J. Chun, 



22 



(2002), |hep-ph/0105157| ; F. Borzumati and J. S. Lee, Phys. Rev. D66, 
115012 (2002), |hep-ph/0207184| ; A. S. Joshipura, R. D. Vaidya and 
S. K. Vempati, Nucl. Phys. B639, 290 (2002), |hep-ph/0203182] . 

[9] M. Drees, S. Pakvasa, X. Tata and T. ter Veldhuis, Phys. Rev. D57, 
5335 (1998), |hep-ph/9712392| ; A. S. Joshipura and S. K. Vempati, 
Phys. Rev. D60, 111303 (1999), |hep-ph/9903435] |; A. S. Joshipura, 
R. D. Vaidya and S. K. Vempati, Phys. Rev. D65, 053018 (2002), 
|hep-ph/0107 204|; V. D. Barger, T. Han, S. Hesselbach and D. Marfa- 
tia, Phys. Lett. B538, 346 (2002), |hep-ph/010826l| . 

[10] S. Davidson and M. Losada, Phys. Rev. D65, 075025 (2002), 
|hep-ph/0010325| ; S. Davidson and M. Losada, JHEP 05, 021 (2000), 
|hep-ph/0005080] |; E. J. Chun, S. K. Kang, C. W. Kim and U. W. Lee, 
Nucl. Phys. B544, 89 (1999), |hep-ph/9807327 . 



[11] H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

[12] G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 (1999), 



hep-ph/9801271|. 

[13] G. F. Giudice, M. A. Luty, H. Murayama and R. Rattazzi, JHEP 12, 
027 (1998), |hep-ph/9810442 . 



[14] L. Randall and R. Sundrum, Nucl. Phys. B557, 79 (1999), 
|hep-th/9810155| . 

[15] KamLAND, K. Eguchi et al, Phys. Rev. Lett. 90, 021802 (2003), 
hep-ex/0212021|; Super-Kamiokande, Y. Fukuda et al, Phys. Rev. Lett. 
81, 1562 (1998), |hep-ex/9807003 ; For a recent review see: S. Pakvasa 



and J. W. F. Valle, hep-ph/0301061 and references therein, 
[16] R. Barbier et al, |hep-ph/9810232; G. Bhattacharyya, hep-ph/9709395. 



23 



[17] Further limits not contained in the reviews |2*1IT6*| can be found in: H. K. 
Dreiner, G. Polesello and M. Thormeier, Phys. Rev. D65, 115006 (2002), 
[hep-ph/Ol 12228 1. B. C. Allanach, A. Dedes and H. K. Dreiner, Phys. 
Rev. D60, 075014 (1999), |hep-ph/9906209| . 

[18] M. Hirsch, H. V. Klapdor-Kleingrothaus and S. G. Kovalenko, 



Phys. Rev. D53, 1329 (1996), hep-ph/9502385|. H. Pas, M. Hirsch 
and H. V. Klapdor-Kleingrothaus, Phys. Lett. B459, 450 (1999), 
|hep-ph/9810" 382|. K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 
75, 2276 (1995), jhep-ph/9506354| . 

[19] S. Rakshit, G. Bhattacharyya and A. Raychaudhuri, Phys. Rev. D59, 
091701 (1999), |hep-ph/9811500| . 

[20] A. G. Akeroyd, M. A. Diaz, J. Ferrandis, M. A. Garcia- Jareno and 
J. W. Valle, Nucl. Phys. B 529 (1998) 3 |arXiv:hep-ph/9707395| . 

[21] A. G. Akeroyd, C. Liu and J. Song, Phys. Rev. D 65 (2002) 015008 
|arXiv:hep-ph/0107218 . 



[22] M. Hirsch, W. Porod, J. C. Romao and J. W. F. Valle, Phys. Rev. D66, 
095006 (2002), |hep-ph/0207334| . 

[23] B. C. Allanach, A. Dedes and H. K. Dreiner, ||arXiv:hep-ph/0309196| 

[24] A. de Gouvea, A. Friedland and H. Murayama, Phys. Rev. D 59 (1999) 
095008 |arXiv:hep-ph/9803481 . 



[25] J. Ellis, S. Rudaz, Phys. Lett. B 128 (1983) 248. J. F. Gunion, 
H. E. Haber, Nucl. Phys. B 272 (1986) 1; Nucl. Phys. B 278 (1986) 
449; [E: Nucl. Phys. B 402 (1993) 567]. 

[26] J. Schechter and J. W. F. Valle, Phys. Rev. D22, 2227 (1980). 

[27] M. Hirsch and J. W. F. Valle, Nucl. Phys. B557, 60 (1999), 



hep-ph/9812463|. 

24 



[28] A. S. Joshipura and M. Nowakowski, Phys. Rev. D 51 (1995) 2421 
|arXiv:hep-ph/9408224 . 



[29] ECFA/DESY LC Physics Working Group, J. A. Aguilar-Saavedra et al. 
|hep-ph/0106315| 



[30] ECFA/DESY SUSY, A. Bartl et al, hep-ph/0301027 



[31] K. Hagiwara et al., Phys. Rev. D66, 010001 (2002) 

[32] T. Behnke, S. Bertolucci, R. D. Heuer and R. Settles, "TESLA: The su- 
perconducting electron positron linear collider with an integrated X-ray 
laser laboratory. Technical design report. Pt. 4: A detector for TESLA," 
DESY-01-011 

R. Brinkmann, K. Flottmann, J. Rossbach, P. Schmuser, N. Walker and 
H. Weise, "TESLA: The superconducting electron positron linear col- 
lider with an integrated X-ray laser laboratory. Technical design report. 
Pt. 2: The accelerator," DESY-01-011 



25 



